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We discuss the force-velocity relations obtained in a two-state crossbridge model for molecular 
motors. They can be calculated analytically in two limiting cases: for a large number and for one 
pair of motors. The effect of the strain-dependent detachment rate on the motor characteristics 
is studied. It can lead to linear, myosin-like, kinesin-like and anomalous curves. In particular, we 
specify the conditions under which oscillatory behavior may be found. 



Understanding the molecular mechanism underlying 
biological motors has recently attracted increasing inter- 
est in biology as well as physics Motor proteins such 
as myosin, kinesin and dynein moving along molecular 
tracks are involved in a wide range of processes essential 
for life, e.g. cell division, muscle contraction, and intra- 
cellular transport of organelles. For many decades ex- 
clusively data from physiological measurements on mus- 
cles H provided experimental information for modeling 
molecular motors |^|. In recent years, a variety of in 
vitro techniques allowed the observation of single motor 
proteins j^] and gave new insights into the basic princi- 
ples underlying their operation. Not only new theoretical 
models for single-molecule motors were inspired by 
these experiments, but also new models for cooperative 
motors gOpQ. 

The theoretical models can follow two different goals. 
Either they are designed to fit as many physiological ex- 
periments as possible by including many (up to six) dif- 
ferent states, or one uses simplified models (mostly with 
two states) in order to extract the generic features of mo- 
tion generation and classify the motors according to their 
properties |f(^[L2|. Latter models fall into two classes, one 
using a specific conformational change (power stroke) in 
the motor molecule ]To|] , the other a ratchet mechanism 
JTH . A striking result of the ratchet models was the pre- 
diction of spontaneous oscillations of cooperative motors 
p3[ , which might explain the oscillatory behavior of mus- 
cles Q. 

Here we discuss the force- velocity relations of a two- 
state model with strain-dependent detachment rates. De- 
pending on the functional form of these rates, the model 
can show a much greater variety of phenomena than pre- 
viously discussed ||[l(|. These include linear, hyperbolic, 
anomalous or kinesin-like force- velocity relations. In the 
two-state model each motor molecule has two long living 
states: attached and detached. This corresponds to the 
model described by Leibler and Huse |Tc| ] when only the 
time limiting steps important for mechanical properties 
are taken into account. Two-state models have also been 
used previously for myosin |3j and kinesin Jt]|| as well as 



in ratchet models |pd| , ^2 15 1. We generalize the two-state 
model by introducing arbitrarily strain-dependent tran- 
sition rates and discrete binding sites. Both extensions 
are crucial for a qualitative and quantitative explanation 
of experiments. 
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FIG. 1. Schematic model for the motor heads running 
along a molecular track. 

The model is defined as follows. Let xh denote the po- 
sition of a (free or bound) head and x m the position of its 
root at the backbone, as sketched in Fig. |l|. Deformations 
of one head can then be described in terms of a harmonic 
potential H — U(xh — x m ~ Xd) with U(x) = ^k m x 2 . 
After attaching to or detaching from the fiber, a confor- 
mational change in the head, described by shifting the 
potential by the distance d m , takes place, i.e. Xd — in 
the attached (A) state and Xd = d m in the detached (D) 
state. This is the first spatial asymmetry in our model 
and constitutes the basic mechanism for the generation 
of directed motion. We assume that the transitions be- 
tween the two states occur stochastically with character- 
istic times t a and td- While the attachment rate t^ 1 can 
be assumed to be constant, there is strong experimental 
evidence fji]|l7j for a strain-dependent detachment rate 
t^ 1 of myosin. 

We will show that different functions t a (x) describing 
strain-dependence of the detachment rate leads to vari- 
ous interesting phenomena, which are the main topic of 
our discussion. In generally t a {x) will be an asymmet- 
ric function, thus bringing a second asymmetry into the 
model. The binding sites are discrete with a separation of 
a = 5.5 nm on actin (8nm on tubulin JTH). Before bind- 
ing to a site, the head position fluctuates due to thermal 
motion. We assume that even in the affine state there 
is a time scale separation between the diffusion and the 



1 



attachment time of the free head. This leads to the prob- 
ability that a head positioned at xq binds to site i given 
by the Boltzmann weight Wi oc exp[— f3U{— xq + ia)] 
p9| . An experimental estimate for the amplitude of 
thermal fluctuations of a free myosin head with data 
from Rcf. Jj"o[ (k m = 0.4pN/nm, d m = 10 nm) gives 
c = \^k B T/k m w 0.3 d m . 

We start our analysis by considering a large group of 
N rigidly coupled independent motors, a situation typical 
for the actin-myosin motor in muscles. Then the fluctu- 
ations resulting from the stochastic operation of single 
motors can be neglected. We set up a Master equation 
for the probability densities $> a (x,t) for a motor being 
attached at x at time t and $d(f) for a motor being de- 
tached. We also need the probability density P(x,x m ) 
that a motor at x m attaches at a distance x from its root. 
This will depend on the actual position of the motor head 
x m with respect to the binding sites [xi = ia). But, since 
generically myosin remains in the detached state for a rel- 
atively long time we may assume the root position before 
attachment to be completely random with respect to the 
discrete binding sites and replace P(x,x m ) by its average 
P(x) = L dx m P{x,x m ) / a. Then, P(x) becomes 
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P(x) = - 



a -f3U(x-d m ) 



a J2 ■ e-0tfO°-dm+3<O 



(1) 



and the probability densities <£> a and $d obey the follow- 
ing Master equations 



(d t -vd x )^ {x ,t) = ^P(x)-^ 



dt* d (t) = 



t d 



dx 



t a (x) 

$a(x,t) 
ta{x) 



(2) 



with normalization $d(t) + J $> a {x,t)dx = N. The force 
produced by the group of motors is given by F(t) = 
J dx$ a (x,t)d x U(x). 
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FIG. 2. a) Examples of a normal an an anomalous F — v 
relation, leading to oscillating behavior, b) A graphical rep- 
resentation of the criterion for the occurance of oscillations. 
If the tangent to the detachment rate as a function of x/d m 
in the point 1 crosses the x-axis right of the point given by 
the duty ratio (d.r.) at zero velocity (as shown above), the 
force- velocity relation is anomalous with a hysteresis around 
v = 0. 

For a constant positive velocity we have to find sta- 
tionary solutions of Eq. H 



G(x,y) 



NjdyG(x,y)P(y) 
vt d + J dx'J dyG(x',y)P(y) 
( f y dx' \ 



(3) 



G{x,y) is a Green's function which can be interpreted 
as the probability that a motor which got bound to the 
fiber at position y still remains bound when its position 
reaches x. 

For a harmonic potential U(x) and a strain- 
independent detachment rate one gets a linear force- 
velocity relation F/N = t a k m (d m — vt a )/(t a + td). ft 
neither depends on temperature nor on the distance be- 
tween the binding sites. More complex functions t a {x) of 
course lead to other forms of the force- velocity relations. 
They can be classified into two groups: the normal ones 
with a monotonously decreasing force for an increasing 
velocity and the anomalous ones, showing hysteretic be- 
havior (Fig. |a). The reason why anomalous relations 
are interesting is that they allow two different velocities 
for the same external force. For instance, in a harmonic 
external potential this leads to spontaneous oscillations if 
the hysteresis spreads over v = 0. Such oscillations were 
first proposed in a two state ratchet model by Julicher 
and Prost Jll],[ll|,?] . Here we show how such a mechanism 
can be implemented in crossbridge model. 

Upon neglecting the discreteness of the binding sites 
and thermal fluctuations, a simple sufficient algebraic 
criterion for the occurance of these oscillations can be 
derived. The zero velocity point certainly lies in a hys- 
teretic range if the slope of the force-velocity relation 
is positive there. Due to the simplification mentioned 
above we set P(y) — S(y — d m ) in Eq. [3] and calculate the 
derivative 



dF 
dv 



N 



kmt a (^d m ) 



v=0 



t a {d Tj 

t d 



td + t a (d m ) 



d m t a (d m )-^ t a 1 (x)\ x=dm ) . (4) 



If its value is positive, the force- velocity relation certainly 
shows anomalous behavior. A graphical representation of 
this criterion is shown in Fig. |[ 

The force-velocity relation as calculated by now de- 
scribes the mean force a group of motors produces when 
moving with a given constant velocity. However, the sit- 
uation is usually reversed and one is interested in the 
mean velocity at a constant force. Of course, both situ- 
ations are equivalent in the limit of large N. But for a 
finite N the motion actually occurs stepwise. This raises 
the question, how the motors remember on which limb 
of the hysteresis they currently move. The quantity that 
actually distinguishes between both limbs is the number 
of currently attached motors. Together with the exter- 
nal force it uniquely defines the velocity. This follows 
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from the fact that the force per attached motor decreases 
monotonously with increasing velocity (Fig. ||), as can be 
seen from Eq. ||. In a finite system the number of motors 
fluctuates and if it passes a threshold value the velocity 
jumps from one stable state into the other. The probabil- 
ity for such jumps is highest if the original state is close 
to the edge of the hysteresis and the number of motors 
low. 
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FIG. 3. An example of the F — v curve for 

t a (x) oc exp(— 2 \x\ /dm) and a low duty ratio. The solid line 
shows the mean force motors would produce at a given ve- 
locity. The dashed line shows the mean number of attached 
motors and the dash-dotted one the mean force per attached 
motor. Seeking for the velocity at a given force one obtains a 
stable (higher v) and an unstable (lower v). A third solution 
is always v — » — oo, but to obtain it one has to include an 
infmitesimally small friction term. 

An example of a function always leading to a normal 
F — v relation is t a (x) = exp(ax). It states that the life- 
time of the attached state is larger for those heads that 
have just gone through the power stroke and produce 
maximum force than for those which have already done 
their work and now pull backwards. As a consequence the 
duty ratio becomes lower at higher velocities. This idea 
has already been used by A.F. Huxley p|. Such a depen- 
dence is needed for an explanation of the approximately 
hyperbolic force-velocity dependence in muscle systems. 
Physiological data by Hill O are perfectly fitted within 
the above analytic results by choosing a d m = 0.55 and 
td ^> t a . Quantized binding sites and thermal fluctua- 
tions are found to play only a minor role (their neglect 
leads to almost the same curve with ad m — 0.58). 

A function that can lead to an anomalous F — v re- 
lation is t a (x) oc exp(— 2 |x| /dm), as shown in Fig. |^. 
For a sufficiently low duty ratio the point v = lies 
within a hysteresis and is unstable if the force is held 
constant. Instead, a positive finite solution is possible 
or a negative with v — > — oo. To obtain the latter from 
our equations, an infmitesimally small friction term has 
to be added. Now if one lets such a group pull against a 
harmonic spring with the other end fixed, the extension 
of the spring oscillates in a sawtooth-like manner with 
flat ascending slopes (motors working against the spring 



force) and very steep descending slopes (the spring force 
pulling motors backwards). Very similar behavior has 
been observed in muscles under some conditions (includ- 
ing low Ca 2+ concentration, which indeed means a low 
duty-ratio) 113] . However, other explanations which sup- 
pose that the oscillations are induced by the regulatory 
system are possible as well pH |. 

So far, our discussion has focused on situations where 
motors are operating in large groups. There is, however, 
a second scenario, where only a few molecular motors 
cooperate at a time, e.g. when kinesin transports vesi- 
cles along microtubules. Modeling them is guided by the 
following experimental observations: A kinesin molecule 
with two heads can move over long distances without de- 
taching from the microtubule p3| . Although it is not yet 
completely clear how the two kinesin heads "walk" along 



the protofilaments |23|, there are good arguments to use 
a model with 8 nm periodicity p4| ] where each binding 
site can be occupied only with one head at a time. Sin- 
gle headed kinesin can move microtubules if cooperating 
in larger groups, but not as fast as double-headed [5s[ . 
The velocity decreases with increasing load almost lin- 
early |^,^{|. Above the stall force kinesin shows back- 
and-forth movement, but does not walk backwards [27|. 
Forward loads can increase the velocity many times [27|. 

From the fluctuation analysis [ p8|j29| it is quite certain 
that both the process of attachment and of detachment 
include an asymmetry between the forward and the back- 
ward direction (a new head attaches in front of the other 
one and the rear head detaches more probably than the 
front one). The low variance |?0| (r w 0.5) is not a priori 
inconsistent with a model where only one of both symme- 
tries is present (e.g. one where steps with 8nm and Onm 
occur with equal probabilities), but then everything else 
in the duty cycle including the dwell times would have to 
be completely deterministic, which does not seem real- 
istic. The behavior at superstall forces j27j additionally 
implies that one of these both asymmetries remains over 
the whole force range, while the other one reverses at 
higher forces. In our discussion we restrict ourselves to a 
model in which the attachment asymmetry gets reversed 
with increasing load while the detachment asymmetry re- 
mains. This approach has already been used by Peskin 
and Oster JtJ] and in a similar way by Duke and Leibler 
||. This, however, does not mean that we consider the 
other case less realistic. 

The central result is again the force-velocity relation. 
Because the velocity is not temporally constant as for 
N — > oo, it has to be calculated directly from tran- 
sition rates for a constant force. Since an attachment 
two sites away from the other head is very improba- 
ble, we take only the attachment rate at the front (f) 
or rear (r) side of the other head into account: R^J r = 
Af exp[±pk m a(d m — F/k m )/2]; Af is chosen such that 
both rates add up to tl . The respective detachment 
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rates are R S J r = t~ x {{±a + F/k m )/2). This gives 



v(F) 



a/2 



'R r d -R{ , 



(5) 



The linear force-velocity curve has led some authors 
p3| , p2| to the conclusion that there is no strain- 
dependence of the detachment rates. This conclusion, 
however, is only valid in a model with continuous bind- 
ing sites |pxj| . Taking into account discrete binding sites 
actually leads to a nonlinear, S-shaped curve if the de- 
tachment rates are strain- independent. Again the sim- 
plest choice is t a (x) — exp(ax). Using a reasonable set 
of parameters the model is able to reproduce the nearly 
linear dependence reported in Refs. [22||{| (Fig. |||) with 
extensions similar to those in Ref. \27\. Fig. |J further 
shows the velocity for a large number of coupled double- 
headed and single- headed kinesin molecules. When com- 
paring them to experiments care has to be taken since the 
pairs are in reality coupled elastically to the backbone, 
which leads to lower velocities. Beside that both curves 
depend very sensitively on the choice of t a (x). Never- 
theless, they show clearly that the "repulsion" between 
heads already causes a significant velocity difference be- 
tween single- and double-headed kinesin. 



iting cases: for a large number of rigidly coupled motors 
and for one pair. In the first case we show how different 
functions describing the strain-dependence of the detach- 
ment rate result in linear, hyperbolic or even anomalous 
force- velocity relations and give a simple algebraic crite- 
rion for the latter. Discrete binding sites play only a mi- 
nor role. For one pair of motors force-velocity-relations 
as measured on kinesin can be reproduced. They de- 
pend crucially on the displacement between the binding 
sites. The model also shows a significant difference be- 
tween single- and double-headed kinesin when operating 
in large grups. 
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FIG. 4. Force-velocity curves for kinesin. Experimental 
data from Ref. at saturating ATP concentration and the 
theoretical curve for a = 0.5 nm -1 , td = 3 ms, t° a — 700 ms, 
dm = 4nm and k m = l.lpN/nm. The middle curve shows 
the result for one double-headed molecule, the upper one for 
many coupled double-headed molecules and the lower one for 
many single-headed molecules. 

Another quantity of interest is the probability for the 
whole molecule detaching from the microtubule during 
one step is Pl = td/t a (F/k m ), yielding 5% at zero 
load and saturating ATP concentration, somewhat higher 
than comparable observations (1.3%) pq]. 

In summary, we have shown that a generalized two- 
state crossbridge model for molecular motors can lead 
to a much larger variety of phenomena than previously 
discussed. We have found analytical results in two lim- 
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